Criticality versus q in the 2 + 1-dimensional Zq clock model 
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Using Monte Carlo simulations we have studied the d = i Zq clock model in two different repre- 
sentations, the phase-representation and the loop/dumbbell-gas (LDG) representation. We find that 
for g > 5 the critical exponents a and v for the specific heat and the correlation length, respectively, 
take on values corresponding to the case q ^ oo, where limq_oo Zq — 3DXY model. Hence in terms 
of critical properties the limiting behaviour is reached already at g = 5. 
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Matter coupled-gauge field theories in 2 -I- 1 dimen- 
sions have come under renewed scrutiny in the context 
of condensed matter physics in the past decade, as effec- 
tive theories of strongly correlated system Concepts 
such as confinement-deconfinement transitions associated 
with the proliferation and recombination of topological 
defects of gauge fields, enter for instance in attempts 
at providing a theoretical foundation for breakdown of 
Fermi-liquid theory in more than one dimension. A 
large variety of such gauge-field theories have been pro- 
posed, and one model of particular interest is the com- 
pact abehan Higgs model[l|,|l,|3)l3- This model consists 
of a compact gauge field coupled minimally to a bosonic 
scalar field with the gauge charge q. In a particular limit 
the dual of this model reduces to a loop-gas representa- 
tion of the global Zq model0,@. This identification has 
been the motivation for the current work, for a detailed 
account of the q dependence of the full theory we refer 
to Ref. HH 

The spin Zq model is a simple planar-spin model, 
where the direction of the spin is parametrized by a 
phase. This phase is restricted to the values 2Trn/q with 
n E Z, and is defined by the following action 
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S — ~p cos ( — (ui 
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The state is specified by the integer variables rii G 
[0, 1,-- - — 1]. Special cases include q — 2 which is 
the Ising model, 9 = 3 which is the three state Potts 
model, and the limit g — > oo which corresponds to the 
XY model. In addition it is easy to see that for (7 = 4 the 
partition function Z(2/3,4) = Z(/3,2) x Z(/3,2). The aim 
of the current paper has been to determine how the crit- 
ical properties interpolate between the well known Ising 
(q = 2) and XY (g 00) limits. We have done this 
by measuring the exponent combination (1 + oi)lv as a 
function of q. 



In d = 2 the model has a quite peculiar phase struc- 
ture, with an intermediate incompletely ordered phase 
(lOP), where the system shows behaviour similar to the 
critical Kosterlitz Thouless phase. Upon further cool- 
ing, the system will order completely into one of the q 
completely ordered states 0, 0- In d = 3 the Zq model 
does not have an lOP, but there are generalisations of 
the model which do [llllioll. 

A related case is that of a globally U{1) symmetric 
theory which is perturbed by a weak crystal field. Using 
RG and duality arguments it has been shown that for q > 
5 the crystal field is an irrelevant perturbation, whereas 
for g < 4 the XY fixed point is rendered unstableTTl. 

It is important to emphasise that we have focused on 
the properties of the Zq model at the critical point. For 
T < Tc the discrete nature of the model will always be 
apparent. A beautiful RG study of the Zq model shows 
how the couplings of the model flow towards a fixed point 
which is ultimately different from the 3DXY fixed point 
in the T ^ limit 

Eq. n is strahtforwardly reformulated as a model of 
an interacting ensemble olinks which either form closed 
loops or originate and termine at point charges. We start 
with the partition function 



Z{/3,q) = ^exp 



(2) 



The first step is to replace the cosine with a quadratic 
potential, this is the Villain approximation . Next, we 
promote the integers nt to real- valued phase variables Oi , 
at the expense of introducing an auxiliary integer field 
Q, which through the Poisson summation formula re- 
stricts the 6i variables to the discrete values allowed by 
original theory. The resulting partition function is then 
given by 
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Zv{P, q) = m j DOY^e^^ - ^ - 27rk)' + iqOQ^ 



{k,Q} 



(3) 



In Eq. 13 {k} is an integer link field living on the links 
of the original lattice, and {Q} is a scalar field living 
on the sites of the same lattice. The prefactor S[/3] and 
effective coupling j3v — Pv{P) must be retained to get 
results which agree with Eq.[21on a quantitative level|l5|. 
however they do not affect the critical properties and 
from now on we will assume (3v = /3, = 1, and omit 
the V index on the partition function. 

In Eq. 121 the Q-field explicitly accounts for the dis- 



crete nature of the Zq model. Setting Q = 0, we re- 
cover the Villain representation of the XY model. Due to 
this similarity, the remaining analysis follows well known 
steps 0, which we briefly include for completeness. A 
Hubbard- Stratonovich decoupling of the quadratic ex- 
pression in Eq. O is performed by introducing an aux- 
iliary field V, thus bringing the partition function onto 
the form 
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Z{P, 1) = J DvD9 exp 



{k,Q} 
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In Eq. ^the {k} summation can be performed, thereby 
restricting the velocity field v to integer values denoted 
by 1. In the term coupling A9 and I, a partial integration 
can be performed, such that only appears in the com- 
bination i6 (Al — qQ), from this we get the constraint 



(Al - qQ) = 0. 



(5) 



At this stage the transformation to a loop gas is com- 
plete, and the partition function is given by 



Z{f3,q) 
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This is a theory consisting of the field {1} living on the 
links of the lattice, and the field {Q} which lives on the 
sites. The field {Q} is sibject the constraint Q = 
0, i.e. overall charge neutrality, whereas the field {1} 
must satisfy the local constraint A • 1 = qQ on all lattice 
points. The latter constraint means that every +Q/ — Q 
pair must be joined by q occupied links, in addition we 
can have {1} excitations which are not nucleated to any 
+Q/ — Q pairs, these must form closed loops. Fig. 
shows a typical configuration for the q — 2 model. 

In the compact Abelian Higgs model considered in 
refs. 1^ the fields {1} and {Q} represent vortices and 
monopoles, i.e. they are the topological excitations of 
the matter-field and gauge-field respectively. That in- 
terpretation does not apply in the current case, but the 



interpretation of the {Q} field is that it maintains the 
discrete properties of the original theory Eq. ^ With 
Q = (the 9 — > oo hmit), Eq. reduces to a loop-gas 
with steric repulsion, this is a well known model with an 
inverted XY transition^3|- Note that the special case 
q = 1 effectively represents no constraint. In this case, 
the theory Eq.Elis noninteracting, and sustains no phase 
transition. For all g > 2 Eq.l^has a phase transition, be- 
tween a phase filled with link segments for (3 > (3c, and 
a vacuum phase which does not contain link excitations. 

We have performed Monte Carlo (MC) simulations of 
the Zq model, using both a phase-representation, Eq. ^ 
as well as the loopgas/dumbbell gas (LDG) representa- 
tion, Eq. |S1 The phase representation is simulated as a 
conventional spin simulation. In the LDG representation, 
the fundamental Monte Carlo moves are represented by 
alternating attempts of inserting a closed loop excitation 
of the 1 field and a dumbbell configuration consisting of 
a +Q/ — Q pair connected with an occupied g- valued 
link (the vertical link to the left in Fig. ^ is an example 
of an elementary dumbbell excitations). For q = 2 the 
(vacuum) excitations of a loop or a +Q/ — Q pair have 
the same energy, while for q > 2 the elementary dumb- 
bell excitations are more expensive than the elementary 
loop excitations, and their relative importance diminishes 
with increasing q. 

The main goal has been to determine how the critical 
properties change with q. The central quantity we have 
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Coupling constant 

FIG. 2: Schematic figure showing third moment of action, and 
how data are extracted for FSS analysis. For further details 
of this method see lal. 
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FIG. 1: A typical LDG configuration for the q = 2 (Ising) 
model. Multiply connected links, like the vertical along the 
left edge have much lower entropy than loop/dumbbell com- 
binations, and hence give a relatively small contribution to 
the partition function. 



considered is the connected third order moment of the 
action jQ] 

ll + O! 



((5-(5)r)oc|/3-/3.r 



(7) 



which recently has been demonstrated to yield surpris- 
ingly good scaling results compared to second moments 
p. When approaching the critical point, the correlation 
length ^ diverges as cx \f} — f3c\~'^ ■ Therefore, in a fi- 
nite system of linear extent L we find that the third order 
moment in Eq. \7\ scales with L as 



{{s-{s)r)^L- 



(8) 



The main advantages of the third order moment in Eq. 
are that (1) good quality scaling is achieved for practical 
system sizes even for models with a < 0, e.g. the 3DXY 
model, and (2) one set of measurements gives both the 
combination (1 -|- a)/^ and —\/v independently P|, al- 
though it is more difficult to achieve high precision on the 
latter. A schematic figure of {{S — {S)y) as a fimction of 
coupHng constant is shown in Fig. |2 and figures 01 and ^ 
show finite-size scaling (FSS) of the peak to peak value. 

We have considered systems of size L x L x L with 
L = 8, 10, 12, 16, 20, 24, 32, 40, 48, and up to 2 • lO'^ sweeps 
over the lattice. In addition to the q — A and q — 5 
presented in figures 01 and ^ we have also studied the 
q values q — 6,8, 12, 16 and 24, ref. @ shows results of 
q = 2 simulations of Eq. jSl We find that the combination 
(1 + a)/v changes abruptly from the Z2 value of 1.763 
[isl l to the XY value of 1.467 [3 when increasing q from 
q ^ 4 to q ~ 5. A further increase of q beyond q — 5 does 
not affect the value oi {1 + a)/^, as shown in Fig. [3 
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FIG. 3: This figure shows the scaling of {{S ~ {S}f} for q = 
4(b) and q = 5{*), the results are obtained using the phase 
representation Eq.|21 The q = 4 results show Z2 scaling with 
{l + a)/u — 1.76±0.05, and the q = 5 curve shows XY scaling 
with {l + a)/u = 1.46 ±0.03. 



That the Zq model is in the XY universality class for 
q > 5 must imply that at the critical point the discrete 
structure is rendered irrelevant for these q values. To in- 
vestigate this point further, we have implemented a sim- 
ple real-space RG procedure, which attempts to probe for 
what values of q the discrete nature of Zq model is rele- 
vant at the critical point. We denote the untransformed 
phases and fields as 6q. The renormalized phase at level 
n -f 1 is given by the block spin construction 



9„+i = atan 



J2k cos On{k) 



(9) 



where the sum over k in Eq. 01 is over the eight spins 
in a 2 X 2 X 2 cube. For q — 2, this transformation is 
clearly trivial, since adding a number of phases and tt 
will still give or tt. However, for q > 2 the effective q* 
will increase with n, and for n —> 00 the resulting block 
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FIG. 4: This figure is similar to Fig. |3 but the resuhs are 
obtained using representation Eq. El The g = 4 results show 
Z2 scaling with (1 + Qf)/z/ = 1.70 ± 0.05, and the q — 5 results 
scale with (1 + a)/^ = 1.47 ± 0.06, i.e. qualitatively similar 
to the results in Fig. El 
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FIG. 6: Histograms of 6 after four rescalings for of the crit- 
ical state. For q — 4, the distribution shows clear signs of a 
discrete background, whereas for q = 5 this is not the case. 
The slow variation in the q — 5 histogram is not commensu- 
rable with a wavelength of 2n/5, and probably only due to 
insufficent sampling. 
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FIG. 5: The exponent-combination {1 + a)/!/ versus q. Note 
how it changes value abruptly as q is increased from q — 4 
to q = 5. The dashed lines are the Ising {Z2) and XY (Zoo) 
values of 1.763 and 1.467, respectively. 



spins can take any direction. 

We next investigate whether the system flows towards 
an infinite value of q* or not under such a RG trans- 
formation. This is tantamount to asking whether the 
discrete structure is rendered irrelevant or not on long 
length scales. To this end, at each iteration step n, we 
have recorded histograms h„ (0) of the phase distributions 
on the lattice, and monitored the manner in which this 
histogram flows under rescaling. By purely visual inspec- 
tion we find that for g = 4 the discrete nature of the Zq 
model persists, whereas for g = 5 it is washed away, this 
is illustrated in Fig. ^ 

To study this RG flow at a quantiative level, we have 
written the phase distribution P„(0„) as a sum of har- 
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an.k cos 
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+ bn,k sm 

9 / \ 1 

(10) 

Here, the coefficient a„^fc in Eq. denotes the fc-th 
Fourier-cosine component at RG level n. Clearly, the co- 
efficient an,q is the interesting component, we have stud- 
ied how this coefficient flows under repeated rescaling. 
For (7 = 4 this coefficent shows critical fixed point be- 
haviour, whereas for q = 5 it flows to zero, even for T well 
below the critical temperature, this is shown in Fig. [7| 

Also the LDG representation Eq.|Bl gives a qualitative 
indication that for q > 5 the discrete nature of the theory 
is irrelevant. In this representation the discrete nature is 
represented solely by the Q excitations, so measurements 
of (IQI) should give a quantitative indication of the the 
importance of the discrete structure. Measurements of 
(IQI) at the critical point give {\Q\) « 0.07, 5.9- 10""^ and 
2.75 -10"^ for q — 2,4: and 5 respectively, whereas the link 
density (|1|) « 0.15 for all q. Hence at g = 5 the discrete 
Q excitations hve been completely frozen out, and the 
tangle is essentially identical to the pure-loop tangle of 
the 3DXY model. 

In summary, we have determined the critical exponent 
combination (1 + a)/!/ in the d = 3 Zq spin model for 
q > 4. Using two different representations we have found 
that for q > 5, the combination (1 + a)/!^ takes a value 
which is consistent with the value taken in the 3DXY 
model. Along with other more qualitative indicators this 
means that at the critical point discrete structure finer 
than g = 5 is irrelavant at the critical point, and the 
long distance properties of the theory are determined by 
the larger symmetry group U{1). These results are in 
accordance with RG studies starting with a U{1) sym- 
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metric theory which is perturbed by a perturbation with 
Zq symmetry. 
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FIG. 7: The flow of the coefiicient for g = 4 and q—'a. 
For q = 4, we see that there is a fixed point at the critical 
point, whereas for q = 5 we see that a„,q flows to zero at the 
critical point. In the figure a„,5 flows to zero also for T <Tc\ 
this is a flnite size effect. This coeffeicient will eventually flow 
to inflnity for sufflciently large systems/low T. 
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